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Abstract. We prove that if O C is bounded and IR^\f2 satisfies suitable structural assump- 
tions (for example it has a countable number of connected components), then W^'^{Q) is dense 
in W^'P{Q) for every 1 < p < 2. The main application of this density result is the study of 
stability under boundary variations for nonlinear Neumann problems of the form 
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1. Introduction 

In this paper we prove a density result for Sobolev spaces defined on two dimensional open 
bounded sets. More precisely, for 1 < p < 2 and C open, bounded and belonging to the 
class ApiM."^) of admissible domains (see Definition 13. l|l . we prove that the Sobolev space W^'^{^1) 
is dense in W^'P{^1). The class Ap{M.'^) contains for example domains whose complements have 
a countable number of connected components or even whose complements are Cantor sets with 
small dimension. 

In the case is sufficiently regular (for example if it satisfies a cone condition), this density 
result is trivial because by means of extension operators and convolutions one can prove that 
C°°(ri) is dense in W^'P{fl). The situation is different when il is irregular: extension operators 
cannot be employed, and the density of C°°{il) in W^'P{il) can fail, as in the case the domain 
contains a crack. Even the density of C°°{fl) in W^'P{il.) proved by Meyers and Serrin j^T] which 
holds for every open bounded set f2 cannot be used because the control on the energy of order 
2 is available only well inside, and can be lost approaching the boundary. In this direction, we 
refer the reader to the paper of O'Farrel |2H| for a counterexample to the density of W^'°°{il) in 
W^'P{n) in the case is too irregular. 
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The main motivation of our density result is the study of stabihty under boundary variations 
for two dimensional nonlinear Neumann problems of the form 



(1.1) 



~diy A{x, Vu) + B{x,u) = inn 
A{x, Vu) • J/ = on (90, 



where A : x g^^^^ ^ . ^2 

X R — > R are Caratheodory functions satisfying standard 
monotonicity and growth conditions of order p (see conditions ()4.3fl - ()4.5|) 'l. Namely we are inter- 
ested in the continuity of the map fl ua, where £ W^'P{n) is the solution of (|l.l|l in (see 
Section^ for the precise sense of the continuity of this mapping). 

The density of W^-^ in W-^'P is a key point to infer stability for problem 1)1. 1|1 from that of the 
linear equation 



(1.2) 



— Au + u = f in n 
#i = on dn. 



Stability results for problem (|1.2|) have been obtained by several authors (see for example |1U) . 
[TT] . [T^ . [T^ . [2j). These results hold in generic dimension N under quite restrictive 

assumptions on and its possible perturbations. For example Chenais proved stability for 
(|1.2|l under a uniform cone condition for the perturbed domains, and this condition excludes 
several interesting cases like those of domains containing cracks which are of interest in fracture 
mechanics. Moreover, the cone condition implies the existence of extension operators, and the 
density of W^'^ in W^'P is trivial, so that the stability of (|1.1|) holds under the same assumptions. 

In dimension N = 2 the situation is different, and restrictions only on the topological nature 
of the domains have been individuated in order to achieve stability for (|1.2(l : this is the reason 
why we are interested in density for Sobolev spaces defined on two dimensional, possibly irregular, 
domains. Bucur and Varchon |Sj consider domains whose complements have a uniformly bounded 
number of connected components and prove that, if fl„ Q in the Hausdorff complementary 
topology (see Section |2] for a definition), we have the stability uq^^ uq if and only if 

meas(il„) — > meas(r2). 

Under strict monotonicity assumptions for A and B, Dal Maso, Ebobisse and Ponsiglione ^Hl 
proved that the same conclusion holds for problem (|l.l|l in the case 1 < p < 2, while for p > 2 
stability is in general false (see Remark 3.7]). The main tool they employ is the Mosco 
convergence of W^^'^'(n„) to W^'P{n) (see Section |21 for a definition) which is equivalent to the 
stability of (|l.l|l for every admissible A and B. The Mosco convergence in the case p = 2 is indeed 
a corollary of the stability result by Bucur and Varchon [Hj. Since they make use of conformal 
mappings, and these are not useful in a nonlinear setting, Dal Maso, Ebobisse and Ponsiglione 
provide a different proof of the Mosco convergence based on nonlinear harmonic conjugates. In 
view of our density result, the Mosco convergence when 1 < p < 2 (and hence the stability result 
for ll.l|l ') can be deduced from the case for p = 2 (see Proposition l4.3|l . 
In Section 21 we consider the nonlinear Neumann problems 



(1.3) 



-divA{x,Vu)+b{x)\u\P^^u = h in fl 
A{x, Vu) -1^ = on dn, 



where b G i°°(R^) is such that 6 > 0, and h satisfies suitable assumptions in order to guarantee 
the existence of a solution. These problems introduce some degeneracy with respect to problems 
(|1.1|) because 6 can vanish on subsets of n with positive measure. As a consequence stability 
cannot be studied in terms of Mosco convergence of suitable functional spaces, because the two 
notions are in general not equivalent (see Remark 5.2]), and so in order to prove stability for 
p. 3(1 . the results of Dal Maso, Ebobisse and Ponsiglione cannot be directly used. 

In the case p = 2, and with A{x,S^) — ^, Bucur and Varchon [Hj proved that if the complement 
of r2„ has a uniformly bounded number of connected components and f7„ — > in the Hausdorff 
complementary topology, then stability holds if and only if 

meas(ri„ n {6 > 0}) meas(fi n {5 > 0}). 
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We prove fProposition l4.4|l that the same resuh holds in the nonhnear case 1 < p < 2. In the case 
p > 2, stabihty does not hold in general (see |^ Remark 3.7]). 

A second application of our density result is to a shape optimization problem, namely the 
optimal cutting of a membrane. The admissible cuts we consider are compact and connected 
sets which contain two given points. The case of a quadratic energy has been treated by Bucur, 
Buttazzo and Varchon in 0. In Proposition 14.51 we prove the existence of an optimal cut for 
a nonlinear energy density /(x,^) with growth of order 1 < p < 2 in ^. Moreover we prove a 
stability result for the associated Euler-Lagrange equation, which is of Neumann-Dirichlet type. 
We remark that in order to establish the existence of the optimal cut and the stability for the 
associated equation, the approximation results of Dal Maso, Ebobisse and Ponsiglione 18 in terms 
of Mosco convergence cannot be used (see Remark l4.7f) . 

Finally, in the Appendix, we show how the arguments of Section |31 provide a new proof of a 
result due to ChamboUe ^] concerning the density of W^-^ in the space LD^'^ of two dimensional 
linearized elasticity. Our approach also covers the nonlinear case LD^'P for 1 < p < 2. 

The main step in the proof of our main result is given by Theorem l3.5l which states the density 
of M^^'^(ri) in M^^'''(17) at the level of the gradients. More precisely we prove that for every 
u G W'^''P{^) we have Vu G H, where 

H ■.= {Vv : ve W^-^{n)} C LP(f],M2)_ 

We use the fact that H = (H^)^, where (■)-'- denotes the orthogonal in the sense of Banach spaces. 
Using Helmholtz Decomposition Theorem, in Lemma 13.41 we characterize H-^ as the family of 
fields ip such that Rip — \7(p with £ W^'P (M^) constant on the connected components of \ft, 
where p' is the conjugate exponent of p and R denotes a rotation of 90 degrees counterclockwise. 
Moreover, using the approximation given in Lemma |3.3I and the fact that fl G ^p(R^), we can 
approximate (j> through functions 0„ G W^'^ (M^) which are constant on a neighborhood of \f2. 
Then the orthogonality of Vu and ip follows by integration by parts. 

The paper is organized as follows. In Section[21we introduce the notation and recall some useful 
notions employed in rest of the paper. Section |21 contains the density result (Theorem 13. 8|l . while 
Section^contains the applications to stability of nonlinear Neumann problems and to the optimal 
cutting of a membrane. In the Appendix we prove the density of W^''^ in the spaces of planar 
elasticity. 

2. Notation and Preliminaries 

In this section we introduce the basic notation and recall some notions employed in the rest of 
the paper. 

If A C R^ is open and I < p < +oo, we denote by LP{A) the usual space of p-summable 
functions on A with norm indicated by || • ||p. W^'P{A) will denote the Sobolev space of functions 
in LP (A) whose gradient in the sense of distributions belongs to LP{A,M.^), and we denote by 
Wq"P{A) the closure in W^'P{A) of smooth functions with compact support. 

If £' C R^, we will denote with meas(i?) its A/'-dimensional Lebesgue measure, and by H°'{E) 
its a-dimensional Hausdorff measure (see Chapter 2] for a definition). Moreover, we denote 
by E'^ the complementary set of E, and by 1^; its characteristic function, i.e., Ie{x) = 1 ii x G E, 
^e{x) =0 otherwise. 

Capacity. Let 1 < p < +oo, and let E C R^. We set 

Cp{E) := inf I J \Wu\P + \u\p dx : u & W^'PIR^), u > 1 a.e. on E^ . 
For the properties of capacity, and its relevance in the theory of Sobolev spaces, we refer the reader 

to uni. 

We say that a property V{x) holds Cp-quasi everywhere (abbreviated Cp-q.e.) on a set £^ C R^ 
if it holds for every x & E except a subset N of E such that Cp{N) = 0. 
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If A C is open, every function u G W^'P{A) admits a quasicontinuous representative, i.e., a 
representative u such that for every e > there exists an open set with Cp{B^) < e and u^axb^ 
is continuous. Throughout the paper, we will identify a Sobolev function with its quasicontinuous 
representative. Notice that for p > N, the continuous representative of u (which exists by Sobolev 
Embedding Theorem) is precisely the quasicontinuous representative. We will use the following 
fact: if ^ w strongly in W^'P(A), we have that up to a subsequence Un —* u Cp-q.e. on A. 

The following lemma will be useful in Section |31 and in the Appendix (a different proof can be 
obtained using the arguments contained in j7| Lemma 5.1, Lemma 5.2]). 

Lemma 2.1. Let u G C(M^), K C connected, and let c G R. If u{x) = c for C2-q.e. x G K, 
then u{x) — c for every x K . 

Proof. By assumption we have that there exists N (- K such that C2{N) — and u{x) = c for 
every x G K \ N . If for every x € N there exists Xn G K \ N such that x„ x, by continuity of 
u we conclude that also u{x) — c and the result follows. 

By contradiction, let us assume that there exists x G N such that x ^ K \ N . Then there exists 
f > such that B{x, r) n (if \ iV) = for r < r. Since C2{N) = 0, by UHl Section 4.7.2, Theorem 
4] we have that 'H°'{N) = for every a > 0, and in particular TC^{N) — 0. As a consequence, 
for every < e < f we can find a covering {B{xi^ri)}i^^i of N such that X^ieN ''i ^- 
B :— UiB{xi, ri) and 

S:={r G]0,f[: dB{x,r) B ^ %}. 

We have that meas(S') < e, so that we can find r < f with dB{x,r) n iV = 0. Moreover, up to 
reducing e, we can assume that \ B{x,r) ^ 0, because otherwise we would get that A^ = [x] 
with x isolated from the rest of A", against its connectedness. 
Let us consider 

ATi ■.^K^B{x,r) and K2 K \ B{x,r). 

Ki and K2 are closed in the relative topology of AT. By construction they are not empty, disjoint 
and such that K — Ki U A'2. But this is against the fact that AT is connected, and the proof is 
concluded. □ 



Hausdorff metric on compact sets and HausdorflF complementary topology. Let A be 

open and bounded in M^. Wc indicate the family of all compact subsets of A by IC{A). IC{A) can 
be endowed with the Hausdorff metric dn defined by 

duiKi, K2) :=max< sup dist(a;, A'2), sup dist(j/,A'i) > 

with the conventions dist(2;, 0) = diam(yl) and sup0 — 0, so that c?/f(0. A") = if A' = and 
dH{9,K) — diam(A) if AT 7^ 0. It turns out that IC{A) endowed with the Hausdorff metric is a 
compact space (see e.g. [20! )• 

In order to treat the stability of Neumann problems under boundary variations (see Section ^J, 
we will use the Hausdorff complementary topology on the family of open sets which is defined as 
follows. Let {fln)neN be a sequence of open sets in M^. We say that ri„ — > in the Hausdorff 
complementary topology if for every closed ball B C we have 

B n il"! ^ B n ri" in the Hausdorff metric. 



The Mosco convergence of Sobolev spaces. In Section 0] we will refer to the notion of 
Mosco convergence of Sobolev spaces in connection with stability results for nonlinear Neumann 
problems. For the reader's convenience, we recall here the definition. 

Let (i7„)„gN be a sequence of uniformly bounded open subsets of R^, and let 1 < p < +00. 
For every u„ G W^'P{fln), let us denote by Unlo„ and by Vu„lf2„ the extension to zero outside 
fin of Un and Vm„ respectively. 
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If is a bounded open set in M.^ , we say that W^'P{fln) converges to W^'P{fl) in the sense of 
Mosco if the foUowing two conditions hold. 

(Ml) Mosco-limsup condition. For every u e W^''P{Vl) there exists u„ G M^^^P(ri„) such that 

Vu„lo„ ^ Vulo strongly in LP(R^,R^) 

and 

Unln„ uln strongly in U'{R^). 

{M2) Mosco-liminf condition. If Uk is a sequence of indices converging to +00, {uk)ki£N is a 
sequence such that Uk e VF^'''(il„^) for every k, and Ukln„^ converges weakly in LP(R^) 
to a function (/?, while Vufcln„ converges weakly in to a function then 

there exists m G W^'P(ri) such that (^a = uIq and $ = VuIq. 

Using a diagonal argument, we have that in order to establish (Ml), it suffices to approximate 
functions belonging to a dense subset of W'^'P{il). This fact will be used several times in Section 
il 

3. The density result 

This section is devoted to the proof of the density of W^''^ into W^'P with 1 < p < 2 on a two- 
dimensional domain which satisfies a suitable structural assumption, for example if its complement 
has a countable number of connected components. Recall that the two-dimensional domain is not 
assumed to be regular (for example it may contain a crack), so that extension operators cannot 
be used. 

First of all, we establish the density result at the level of the gradients (Theorem 13.511 . The 
extension to the full result on Sobolev spaces fTheorem lS.SI) is then obtained through a truncation 
argument. 

The class of admissible domains we consider is given in the following definition. 

Definition 3.1. (The class ^p(M^) of admissible domains) Let 1 < p < 2, and let C 

be open and bounded. Let {Ki}i^j be the family of the connected components ofQ"^. We say that 
n belongs to the class Ap{R'^) of admissible domains if for every i £ I there exists Xi G Ki such 
that setting E := {xi, j G /} we have 

(3.1) n^^P{E) = 0. 

Notice that the class 0;(K^) of two-dimensional domains such that their complements have at 
most I connected components (which is relevant for stability of nonlinear Neumann problems, see 
Section r4.1|l is contained in ^p(M^). Moreover Ap{M.'^) contains domains such that fl'^ has a 
countable number of connected components, or even an uncountable number provided that there 
exists a suitable selection E of {Ki}i(zi with zero Hausdorff measure of order 2 — p. We remark 
that condition H3.1fl is not referred to the connected components Ki of Jl'^ but to a selection E of 
{Ki}i^j: in particular it can be meas(iiri) > (not only for the unbounded connected component). 

The following lemmas will be useful in the proof of Theorem 13.51 

Lemma 3.2. Let A CR"^ be open, and let u G W'^''^{A) with q > 2. Then we have meas(u(_B)) = 
for every E C_ A such that Tii-^ (E) — (in the case q — +00 we mean Ti. (E) — 0). 

Proof, li q — +00, the result follows because u is a locally Lipschitz function and meas(/(C)) = 
TO-{f{C)) < LH^{C) for every L-Lipschitz function / and every set C (see ^12, Theorem 1, Section 
2.4.1]). 

In the case 2 < g < -l-oo, we follow the approach that Marcus and Mizel [23 developed to 
deal with iV-property of Sobolev transformations (see [201 for a description of the problem of 
TV-property, and [23 Theorem 5.28]). 

By Sobolev Embedding Theorem u is a Holder continuous function. Moreover, for any square 
Qr A oi side r we have 

(3.2) \u{x) -uqJ< C,||VM|U,(Q^,R2)ri-2/9, 
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where uq^ denotes the average of u on Q^, and Cq depends only on q. From H3.2|l we deduce that 
u{Qr) is contained in an interval Iq^ of length at most 

Iq^ :=2C,||Vu|U,(Q^,R2)ri-2/'?. 

Ijet E CA be such that (E) = 0, and let us fix £ > and 5 > 0. Since n^{E) = 0, we 
can find a covering J- = {Qriixi)}ieN of E with Qnixi) C A, 

(3.3) Y.^r<e 

i=0 

and such that 

By Besicovich Covering Theorem (see ^1 Section f.5.2, Theorem 2]) there exist m families 
J-i,J-2,-.-,Tj,.-.j J-m C of disjoint squares {Qn j such that 

EC\J UQ...,(x„). 

j=l i=0 

By Holder inequality and by 1)3. 3|) we deduce that 

j=l 1=0 j = l 1=0 



1-2/9 



(rn +c>o \ / rn +00 _, 

EEii^-iii^(Q.,(...),«^) EE-ff 
J = l i=0 ' / \j=l ^=0 

q-l 

< 2Cqm||Vu||i,(A^R2)£ 1 

so that 

Hj(7/(i;)) < 2C,m||VM|U,(^^R2)£^, 

where 7{g{E) denotes the (1, (5)-Hausdorff pre-measure. Since H^{u{E)) = lim^^o ^K^C-^))' 
TC'^{u{E)) = meas(w(i?)), we conclude that 

q-l 

Taeas{u{E)) < 2Cgm|| Vu||L9(^^R2)e i . 
Since e is arbitrary, we deduce that meas(M(i?)) = 0. □ 

Lemma 3.3. Let cj) S W^^p{R^) D C"^(R^) with p e [1, +oo]. Let K CR^ be such that (j){K) is 
compact and meas((/)(i^)) = 0. Then there exists 0„ e VV^'^iR^) H C°(R^) with 

(3.4) (j)n^(j) strongly in W^'P{M.^) if I < p < +oo, 

(3.5) (-/-n, V0„) ^ (0, V</>) weakly* in L°°{R^,R^+^) if p ^ +oo, 

and such that 0„ is locally constant on a neighborhood of K, i.e., V0„ = a.e. on a neighborhood 
An ofK. 

Proof. By assumption C := (j){K) is compact and such that meas(C) = 0. Let us set 

Cn:=\yeR: dist(2/,C)<- 

and 

ry 



Tn{y) / '^R\cjs)ds. 
Jo 

Since meas(C„) — s- as n — > +oo, we get that 

(3.6) Tn Id pointwise. 
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Moreover, T„ is 1-Lipschitz, T^^ — a.e. on C„, and 

(3.7) T,;(?/) 1 for a.e. j/ e R. 
Let us set 

(3.8) (t>n ■■= Tn O (j). 

We have that (j)n S W^''p{M.^) n C°(K^), and by the Chain Rule Formula for Sobolev functions 
(see for instance |31 Theorem 3.99]) we get for a.e. x e 

(3.9) V0„(x) =T^(0(a;))V0(x) 

(recall that V(/) = a.e. on 0^^(C) since C has zero measure Proposition 3.92])). 

In view of (|3.9I) . V^n = on yl„ 0^^(C„) which is a neighborhood of K. Moreover, by (|3.6|l 
and H3.7(l . we have that (|3.8|l and (|3.9() imply that 

0„ — > and V(/)„ — » V(/) a.e. on M^. 

Since \4>n\ < \<P\ and |V0„| < \V(f>\, we deduce that (|?^ and hold. □ 

The following lemma is very close in spirit to jl8l Lemma 3.6] 

Lemma 3.4. Let C be open and bounded, and let q>2. Let ip G L'(ri,R^) be such that 

-0 • Vudx = for every u G W^''^{n). 



Then there exists (j) G Vl^^'''(R^) constant on the connected components of Q,'^ (in the case q = 2 
constant C2-quasi everywhere) and such that 

where R(a,b) := (— &, a) denotes a rotation of 90 degrees counterclockwise. 

Proof. Let us denote by Ki, i E I, the connected components of il"^, and let Kq be the unbounded 
one. 

Since e L^(ri,R^), by Helmholtz decomposition of the space L^{fl,M.'^) (see |23 Theorem 
1.1, Chapter III]), there exists ipn G C^(J7,R^) with div^n = and 

(3.10) V-n^V- strongly in i2(f7^R2)^ 

By setting ipn = outside il, we can consider ipn as defined on the entire R^. Let us consider 
ipn '■— Ripn- Since R^ is simply connected, and (pn has zero-curl, we have that there exists 
(f>n e C'°°(R2) such that 

Vcpn = Ripn- 

In particular V(pn = on a neighborhood A" of ^l'^, so that for every i G / there exists c" G R 
such that 

(3.11) </)„ = c" on a neighborhood A" of Ki. 

Since (pn is well defined up to a constant, we can assume that (pn = on Kq. Let D be a disk 
centered at the origin and such that fl C D. By H3.1()(l we deduce that there exists (p G Wq''^{D) 
such that 

(pn (p Strongly in Wq''^{D). 

We have that 

V(p:^ Rip e L'i{n,R^). 

We deduce that <p G Wq^'*(£'), and in particular (p G W^'''{M.'^). Since up to a subsequence </>„ (p 
C2-q.e., from H3.11|l we deduce that there exists q G R, i G /, such that 

(3.12) (p = Ci C2-q.e. on Ki. 

In the case q > 2, we have that (p is Holder continuous by Sobolev Embedding Theorem. So by 
Lemma [2 .11 we get that (I3.12|) implies that (p is constant on Ki, and the proof is concluded. □ 

The following theorem contains the density result for the gradients. 
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Theorem 3.5. Let 1 < p < 2, and let fl G Ap{M.'^) be an admissible domain. Then for every 
u G W^'^{Q) there exists {un)n&i sequence in W^''^{VL) such that 

strongly in LP(fi,R2). 
Proof. Let Ki, i € I, he the connected components of il"^. Let us consider 

(3.13) H ■.= {Vv : ve W^'^{n)} C LP{n,R^). 

In order to prove the density result, it suffices to check that for every u e W^'P{fl) we have 

Vii e i?, 

where H denotes the cfosure of H in the L^-norni. Since 

where (•)^ denotes the orthogonal space in the sense of Banach spaces, we have to check that 

(3.14) VueiH-^)-^. 

Our strategy to prove (|3.14|l is the following. Firstly we characterize the functions V' G H^, 
and then we prove that for every u G W^'P{n) we have the orthogonality condition 

(3.15) / ^p-Vudx^O. 

Jo, 

Step 1: Characterization of H^. Let -0 G {^^^^)^ where p' > 2 is the conjugate 

exponent of p (p' — if p g]1,2[, p' = +oo ii p = 1). By definition of we have that for 
every v G W^'^(ri) 

/ -0 ■ Vt; = 0. 
Jn 

By Lemma IX^ we deduce that there exists 4> G W^'P' (M?) with V0 = R^Jj {R{a,b) := {-b,a) is 
the rotation of 90 degrees counterclockwise), and such that for every z G / 

(3.16) (f> — Ci on Ki 
for suitable Ci G M. 

Step 2: Checking the orthogonality condition. In order to conclude the proof, it suffices to 
check that (|3.15|l holds for every u G W^'P{fl). By Step 1, we need to check that 



(3.17) / ip-Vudx^- RVip -Vudx = Q, 

Jq Jn 

where G W^'P'(R^) satisfies for some q G M, « G /. 

Notice that is compact. Moreover, since G ^p(R^), there exists a selection E of the 

connected components Ki of fl'^ such that 

n^^{E) ^n^-p{E) ^0 ifi<p<2 

and 

n^{E)=0 ifp=l. 

By (|3.1fi(l we get = 4>{E), and by Lemma Hi 21 we have that meas(0(rj'^)) = meaa{(j){E)) = 0. 

Applying Lemma [3.31 there exists 0„ G W^'P (M^) such that 

(3.18) strongly in Ty^'P'(M2) if 1< p < 2, 

(3.19) (0„, V0„) ^ (0, V0) weakly* in L^{R^,R^) if p = 1, 
and 

(3.20) ^(pn = on a neighborhood An of fi'^. 
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Notice that RV4>n is divergence-free. Up to reducing A.^, we can assume that 57 \ A.^ is regular, 
and that the support of R\7(f)„ is contained in \ An- Then we have 



R\7(j) ■ Vu dx — lim / i?V(/>„ • Vm dx — lim / RVipn ■ Vw dx, 

and integrating by parts we deduce that 

lim / i?V0„ • Vu da; = lim / div {RS/(f>n) udx — 0, 

SO that (|;-{.17|l is proved, and the proof is concluded. □ 

Remark 3.6. As mentioned in the Introduction, the density result given by Theorem 13.51 fand 
the similar result for Sobolev spaces Theorem I3.8|l is useful to establish a link between stability 
results for linear and nonlinear Neumann problems. Since stability results usually hold under the 
assumption of a uniform bound on the number of the connected components of the complements 
of the varying domains (see Section ing , the case O'^ has a finite number of connected components 
is the relevant one for the applications. 

In this case, the existence of the function 0„ satisfying conditions H3.18|l and (|3.2()|l in Step 2 
of the proof of Theorem 13.51 can be established more directly without using the approximation 
Lemma 13.31 as follows (the case p = 1 is usually not considered in the study of nonlinear Neumann 
problems in view of a lack of compactness of W^'^). 

Let Kq, Ki, . . . , Km be the connected components of rj"^, where Kq is the unbounded one. Let 
us consider G C°°(]R'^) and e Cj?°(R^), i = 1, . . . , m such that = 1 on a neighborhood of 
Kq, = 1 on a neighborhood of Ki, and 

supp(^/i) n supp(Cfe) = iov h^k. 

By Theorem 9.1.3] for every i = 0, 1, . . . , m we can find < G C°°(R2) with 

0^ = Ci on a neighborhood of Ki 

and such that 
Setting 



strongly in ^^^'^'(M^). 




we get that (|XTH|) and hold. 

Remark 3.7. In the proof of Theorem 13.51 we used the assumption that belongs to the class 
^p(R^) in order to apply the approximation Lemma f3. 31 and recover the functions satisfying 
(|rTH|) . l|TT^ and (ESni- Lemma requires that meas((/)(f]'=)) = 0, and for Q. € ApiR"^) every 
function (p S W^'^ (M^) constant on the connected components of ^l'^ is such that meas((/)(ri^)) = 0. 
In particular this is the case for the functions we need to approximate, that is (j> G VF^'^ (R^) such 
that i?V0 e H^. 

We do not know if for a general ft we can have meas((/)(ri'^)) = for any </> S W^ ^ (R^) 
determining an clement of H-^. For such a (j>, by Step 1 (and in view of the proof of Lemma l3. 4(1 . 
we have that there exists a sequence of smooth functions 0„ such that 

(3.21) V0„ = on a neighborhood of il'^ 
and 

(3.22) ^ strongly in W^'^{R^). 

In particular (j)n{^'^) is finite so that meas((/)„(ri'^)) = 0. If this always implies that in the limit 
meas(0(r2'^)) = 0, the fact that (p is energetically more regular than i.e., (f> € W^ ^ (R^), plays 
an essential role. 

We can consider indeed the following example which shows a sequence (0n)neN of smooth 
functions satisfying and but with (j) G W^^'^{R^) H C{R^) and such that (/)(f7'=) is the 



10 



A. GIACOMINI AND P. TREBESCHI 



interval [—1, 1]. Moreover ^l'^ can be chosen such that its connected components admit a selection 
E with dimension zero, i.e., H"{E) = for every a > 0, so that 17 G Ap{M.'^). This example heavily 
relies on a construction proposed by Maly and Martio j23l in connection with the A^-property of 
Sobolev transformations. 

Let us consider the square Q :=] - 2, 2[x] - 2, 2[ in R"^, J {(t, 0) : -1 < t < 1}, and a„ \ 0. 
Since a point has C2-capacity zero, there are functions Um G C°°(R^) such that 

u,n ^ strongly in W^'^{M.^) 

and such that < < 1, £ int{M„i = 1}, and u„i — outside the ball i?(0, 1). Let zi, Z2 G J 
and ro > be such that the balls _B(zi,ro) and _B(z2,ro) are disjoint. Let us set 

1 f X - Zi\ 1 



V ''0 / 2 V ro 

The functions € C°°(R^) are constructed as follows. Let be the constant function equal 
to 0. If n > 1, let us divide the interval / [^l, 1] in n intervals /" of length 2/n; we can find 
points xf £ J and a radius r„ so small that nr"" — s- and (t>n-i maps B(xf,rn) to the middle 
point of /f . Let B„ ULi ^(^^r.^n), 

and let m„ be such that 
We set 



2 "^^5m (^"~7~^j if N ~ I < r„ for some i 
otherwise, 

Il^m„,n||wi.2(][j2) < 2 



[>n ■= (Pn-1 



and we denote by (j) the strong limit in W ' (S. ) of {(f>n)neN, which is by construction a Cauchy 
sequence. Notice that cj) S VF-^'^(R^) n C(M^), and that the convergence is also uniform. 

Let := Q \ f^nm^"-- have that Q'^ = Q'^ L) HnGN-^"- Since nr"" 0, we have that 
7i"(n„gN^n) = for every a > 0. As a consequence, the connected components of V,'^ admit a 
selection E such that Ti^iE) = for every a > 0. In particular G ^p(R^). 

By construction we have that </)„ is constant on a neighborhood of Q'^ but, since 0„ — > </> 
uniformly, it is easy to see that 4>{^'^) = [—1,1]. Clearly (p cannot belong to W'^''^{M.'^) for some 
q > 2, because otherwise its Holder continuity would imply meas(<^(ri'^)) = 0. 

We are now in a position to prove the main density result of this paper. 

Theorem 3.8. Let 1 < p < 2, and let G Ap{M.'^) be an admissible domain. Then W^'^{^1) is 
dense in W^^P{n). 

Proof. The main ingredients in the proof are Theorem 13.51 and a truncation argument. It is not 
restrictive to assume that is connected, because we can work on each connected component. 

Let u G W^'P{Q). The density result will be proved if we show that for every £ > we can find 
{un)n<£N sequence in W^ '^{n) such that 

(3.23) limsup ||u„ - u||H/i,P(a) < e^, 

n — >+oo 

where ^ as e ^ 0. 

It is not restrictive to assume that 

uew^'P{n)nL°^{n). 

In fact, if fc > and 

Tk(u) min{max{u, — fc}, fc}, 
we have Tk{u) G W^-P{n) n and Tk{u) u strongly in W^^P{n) as fc ^ +oo. Then if 

(|3.23|) holds for Tk{u), by a diagonal argument it also holds for u. 
Let A CC ^l, A regular and connected, such that 

(3.24) \\Vu\\l,^^^^^ + \\u\an\A\<e. 



A DENSITY RESULT FOR SOBOLEV SPACES IN DIMENSION TWO 



11 



By Theorem 13. 51 there exists G W^'^{il) such that 

(3.25) Vvn^yu strongly in LP (17, R2). 
We claim that, up to adding a constant to u„, we can assume that 

(3.26) Vn^u strongly in W^-p{A). 
Let us set 

(3.27) Un ■■= min{max{w„, -||u||oo}, ||u||oo}- 
Notice that u„ G W^'^{n), 

(3.28) Un^u strongly in Tyi'P(^), 
and that 

(3.29) |Vu„| < \Vvn\ a.e. in ft. 

In view of l|I?7jl . and (IX^ . we deduce that 

limsup ||u„ - u||^i,p(f^, < limsup ||u„ - u||^i,p(^) + limsup ||u„ - uf ^ 



limsup||M„-ur , ^ <limsup2P M / 



_ yUn\P +\yu\P +\U„\P +\u\Pdx] 
Q\A J 



so that (|3.23|) is proved. 

In order to complete the proof, let us check that claim H3.26|l holds. If 



A 



1 

since A is regular, by Poincare inequality we have 

Vn — Vn ~ c„ is bouudcd in W^'P{A). 

Moreover, by the compact embedding of W^'P{A) in LP(A), there exists v e L'''{A) such that up 
to a subsequence 

Vn ^ V strongly in L''{A). 

Since Vw„ = Vv„ on A, and in view of (|3.25|) . we get that Vw = Vu in the sense of distributions 
on A. We deduce that v £ W^'P{A), 

Vn ^ V strongly in W^'P{A), 

and since A is connected 

V = U + CA 

for some constant ca £ K. If we set Vn '■— Wn — c„ — ca we get 

Vn u strongly in W^'P{A), 
so that claim (|3.25|l is proved. □ 

Remark 3.9. (The case unbounded) The density of W'^''^{n) into W^'P{n) when 1 < p < 2 
holds also in the case is unbounded but there exists r„ — > +oo with Qn '■= ilr\B{0, r„) G ^p(R^). 
In fact, if u e W^'P{n), we have u £ W^^PiQn) so that there exists e W^^^i^n) with ^ u 
strongly in VF^'P(S1„) as fc +oo. If Xn is C°° function with < Xn < 1, Xn = 1 on 5(0, r„/2) 
and Xn = outside B{0, r„), in order to conclude it suffices to choose m„ :— Xn G VF^'^(r2) for 
kn sufficiently large. 
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4. Applications 

In this section, we give some applications of TlLeorein l3.8l to stability under boundary variations 
of nonlinear Neumann problems, and to the optimal cutting of a membrane. 

Since we work on domains which are not assumed to be regular (for example they may contain 
cracks), we will use Deny-Lions spaces. They are defined as follows. Let f2 be an open subset of 
R^, p e [1, oo[, and b S L°°{R^) with 6 > 0. Let us set 

Cl'P{n) := {u e iL(^) : e LP{n,R^), [ \u\Pbdx < +00} . 

Jn 

We say that uTZbV if 

/ [|Vu - Vv\P + b\u-v\P]dx^ 0, 
Jn 

and we set 

(4.1) Ll^P{n):=£l^P{n)/TZ, 

endowed with the norm ||w|| :— || VwH^pj-fj kn-j + (J^ \u\Pbdx)^^P. L^'P{fl) is the Deny-Lions space 
of exponent p and weight b. In the case 6 = 0, it is usually denoted by L-^'P{il). 

Notice that W^'P{n) C Ll'^in). In the case 6 > c > and is Lipschitz, we have that equality 
holds, while if b vanishes on subsets with positive measure or $7 is irregular, the inclusion can be 
strict (see for example [221 Section 2.7]). Moreover W^'P^fl) is always dense in Ll'P{n), as one 
can check by truncation. As a consequence, in view of Theorem 13. 81 we have the following density 
result. 

Proposition 4.1. Let G Ap{M.'^) be an admissible domain (see Definition and let b G 
L°°(]R2) such that b>0. Then W^^'^{n) is dense m Ll'^in) forl<p<2. 

Let 

(4.2) Oi(M^) {A C open : \ A has at most / connected components}. 

For every u G L^'^(f2), we denote by VuIq and uln the extension to zero outside ft of Vm and 
u respectively. We will use the following proposition due to Bucur and Varchon (see [S] and |H1 
Theorem 4.1, Remark 5.2]), which is a sort of Mosco limsup condition (see Sectional) for the 
spaces L^'^. 

Proposition 4.2. Let r2„ be a sequence in Oi (R^) converging to fl in the Hausdorff complementary 
topology (see Section\^ and such that 

meas(fi„ Ci {b > 0}) meas(f2 n {6 > 0}). 

Then for every u G i^'^(51) there exists u„ G L^'^(f2„) such that 

Vm„1o„ Vulo strongly in ^^(R^R^) 

and 

Wnlf2„ uln strongly in L^(R^), 
where i^(R^) denotes the L'^ -space on with weight b. 

Propositions 14 . l1 and 14 . 2l will be our main tools in dealing with stability of nonlinear Neumann 
problems and with the optimal cutting of a membrane. 

4.1. Stability of nonlinear Neumann problems under boundary variations. Let p G 

]1, +00 [, and let A : R^ x R^ ^ R^ and B : x M ^ R be two Caratheodory functions such that 
the following conditions hold: there exist a G L^' {M?) {p' := ^), P G i^(R^), < ci < C2 such 
that for almost every cc G R^ and for every ^, ^1, ^2 G R^ with ^1 7^ ^2 

(4.3) (^(x,a)-^(2;,6))(a-6)>o, 



(4.4) 



\A{x,0\<<^{x)+C2\£,\ 
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(4.5) A{x,0-i>m+(^i\i\''- 

We assume that B satisfies (|4.3ll . (I4.4|l and l|4.5|l for almost every x £ M^, and for all ^, ^2 G 
with ^1 ^ ^2- 

Let J7 be a bounded open subset of R^. We are interested in the stability under boundary 
variations of f2 of the elliptic equation 



(4.6) 



-div A{x,Vu) + B{x,u) ^ in il, 
A{x, V-u) • 1/ = ondn, 



where v denotes the outer normal to dfl. Since we do not assume any regularity on A, B and on 
the boundary of fi, we intend (|4.6() in the usual weak sense of Sobolev spaces. More precisely by 
a solution of problem H4.6|l we mean a function uq G W^'P{n) such that for every test function 
ifi e W^'P{n) we have 



/ [A{x, \7un)\7(p + B{x, un)^] dx = 0. 



Existence and uniqueness of a solution to 1)4.6(1 follow by well known results on nonlinear elliptic 
equations with strictly monotone operators (see for instance \Z2\). 

Let (r2„)„gN be a sequence of uniformly bounded open sets in . Wc say that Q. is stable for 
the Neumann problems 1(4. 6|l along the sequence (f2„)„gN if 

wn„ln„ ^ unln strongly in L^iB?) 

and 

Vuo„lo„ ^ Vuolo strongly in (R^M2). 

Dal Maso, Ebobisse and Ponsiglione |18l Theorem 2.3] proved that Vl is stable for problem ((4.6|l 
along (r2„)„gN for every admissible A and B if and only if the space W^''''{Q,n) converges in the 
sense of Mosco to W^'P{VI) (see Section |21 for a definition). 

If ri„ is in the class ©((R^) defined in (|4.2|) . Bucur and Varchon 8 proved that if f2„ 57 in 
the Hausdorff complementary topology (see Section [3 for a definition) then VF^'^(f2„) converges 
in the sense of Mosco to M^^^^(f2) if and only if meas(0„) — >• meas(J7). 

Dal Maso, Ebobisse and Ponsiglione jTH] extend this result to the case 1 < p < 2 using a 
technique of nonlinear harmonic conjugates, and they prove that in the case p > 2 the result is in 
general false. 

In the following proposition, using Theorem 13. 81 we prove that the result of JS] can be deduced 
directly by that of [Hj. 

Proposition 4.3. Let (r2„)„gN be a sequence of uniformly hounded sets in Oi{M.'^) such that r2„ 
converges to il in the Hausdorff complementary topology, and let 1 < p < 2. Then W^'P{iln) 
converges in the sense of Mosco to W^'P{n) (and hence problems (|4.6f) are stable) if and only if 

(4.7) meas(r2„) meas(r2). 

Proof. Let us assume that the Mosco convergence holds. Let ^ e R^ with |f | — 1. Let us consider 
u e W'^'P{il) such that u{x) :— ■ x. By (Afl)-condition, there exists ii„ e W^'P{iln) such that 

Vu„lji„ -> Vwln strongly in LP{m.'^,R^). 

Since |Vu„lf2^ — Vwljij = 1 a.e. on f2 \ i7„, we get 

(4.8) limsupmeas(r2 \ r2„) < lim / iS/unln^ — Vul^]^ dx — 0. 

Let us consider m„ G W^'P{fln) such that u„(x) = ^ ■ x. By (M2)-condition, up to a subsequence 
we have that there exists u G W^'P{ft) such that 

Vm„1o„ ^ Vulo weakly in LP{R^, M^). 
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Then, if _D is a disk containing ri„ for every n we get 

(4.9) ^ lim nieas(ri„ \ fi) = lim / Vunln^^DXn dx ^ / \7ulnlD\ndx = 0. 

Combining and we get that holds. 

On the contrary, let us assume that (|4.7|l holds, and let us prove the Mosco convergence of 
W^'Pifln) to w^'P{n). 

Concerning condition (A/2), let Uk G W^'P{^lnk) be such that 

Vufclo„, weakly in LP {R^,R^) 

and 

Wfcla^j, ^ </? weakly in i''(R^) 

for some $ G LP(R^,K^) and ip £ LP(R^). Clearly, since ri„ converges to in the Hausdorff 
complementary topology, we have that <& = Vip on fi, so that u :— {(p)^^ G W^^'P(ri). In order to 
conclude that (M2) holds, we have to prove that $ = VuIq and (f = uIq. Since f7„j. ^ in the 
Hausdorff complementary topology and meas(f2„) meas(r2), we have 

lo„, ^Ifi strongly in ii(R2). 

Then for all rj e C;?°(R2,M2) deduce 

/ ^-rjdx^ lim / [Vufclfj^ ] • ry = lim / [Vufclo„. ] • [f^ln^, ] da; 

$-?7lndx= / ^ln-r]dx^ / Vuln-r/dx 



so that $ = VuIq. Similarly we can prove that ip ~ uIq. 

Let us prove condition (Ml). Since it is sufRcient to approximate functions in a dense subset of 
W'^-P{n), since ©/(M^) C Ap{R'^), we can consider in view of Theorem IXFl functions u e W^-'^{n). 
Then by Proposition l4.2l fwith 6=1) there exists u„ € W^''^{nn) such that 

||Unlf2„ - m1o||l2(r2) + ||Vu„lo„ - Vm1o||l2(r2^r2) -> 0. 

Since W^''^{^ln) Q W^'P{nn), and since the L^-norm is stronger that i^'-norm on bounded domains 
(1 < p < 2), we deduce that (Ml) holds, and the proof is concluded. □ 



Let us now consider the following nonlinear Neumann problem 
(4.10) 



-div A{x,\7u) +b{x)\u\P^^u ^ h in fi, 
A{x,\7u)-v = ondfl, 



where A is a Caratheodory function satisfying conditions H4.3|l . 14. 4|) . (|4.5|l . A{x,0) — for a.e. 
X £ R^, & £ L°°(R^) and b > 0. In order to guarantee the solvability of (|4.1U|) . we assume moreover 
that h = bf + g with f,g^LP (R^), supp C fj, and g dx = for every connected component 
C of fi. We are interested in problem 14.10|l because it introduces some degeneracy with respect 
to problem 14.6|l as b can vanish on regions of fl with positive measure. 

By a solution of l|4.10|l we mean a function uji £ Lf^'P{Q) (or more precisely an equivalence 
class, see (|4.1f) for a definition) such that for every ip £ Ll'P{n) 

[A{x,'Vun) -Vip + b{x)\un\''~'^unp]dx = / hipdx. 
n Jn 

Notice that the integrals appearing in the weak formulation of (|4.1U|) are well defined: in particular 
notice that, if [/ is a regular open set such that suppg C C/ C J7 C fi, and (p denotes the average 
of iy9 on J7, by Holder and Poincare inequalities we get 



gp dx 



n 



gpdx 



u 



g{ip - (p) dx 



u 



< 



C|l3llLP'(R2)l|V(y5||LP(0)- 
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The existence of a solution of (|4.1(J|I can be estabhshed minimizing on L^'^(il) the functional 

F{u) := / [A{x,Vu) ■yu + b{x)\u\P -hu\dx 
Jn 

by means of the Direct Method of the Calculus of Variations. Uniqueness of the solution follows 
by strict convexity of F. 

We say that is stable for the Neumann problems H4.10|l along the sequence (r2„)„gN if 

lim / |Vuo„lj2„ - Vunlfil^ + 6(x)|-un„lo„ - ttfilfj^ rfa: = 0. 

n — ^+00 / 

The stability of Neumann problems (|4.10|l has been investigated by Bucur and Varchon in ^ 
in the case p — 2 and A{x,S,) = ^ (but it easily generalizes to = with a giving the 

correct coercivity). 

The main interest in the stability of H4.10|l is that, since b is not assumed to be strictly positive, 
stability is not equivalent to the Mosco convergence of i^' (a^) to Ll'in) (see Remark 5.2]). 
As a consequence, passing to the nonlinear setting with 1 < p < 2, an approach to stability in the 
line of Dal Maso, Ebobisse and Ponsiglione ^Hl based on Mosco convergence cannot be directly 
used in this situation. 

Let (J7„)„gN be a sequence of uniformly bounded open sets in M^. In the case p — 2, Bucur and 
Varchon ^ proved that, if Sl„ S Oi(R.'^) and ri„ — *■ in the Hausdorff complementary topology, 
then stability of (|4.1U|) holds if and only if meas(f7„ n {6 > 0}) meas(57 n {6 > 0}). Proposition 
14. II permits to extend this result to problems (|4.10() . 

Proposition 4.4. Let {Qn)neN be a sequence of uniformly bounded open sets in Oi{R'^) converging 
to D, in the Hausdorff complementary topology. Then is stable along (i7„)„gN for the Neumann 
problems (|4.10|l if and only if 

(4.11) meas(17„ n {6 > 0}) -> meas(fl n {6 > 0}). 

Proof. Let us assume that stability holds. Then if we choose / = 1 and g = so that h — b, we 
deduce that 

uo„ = ln„ and uo = lo 

and that 

(4.12) ln„ lo strongly in Ll{R^), 
where L^(R^) denotes the L^-space with weight 6. In particular 

(4.13) meas([rJ„Ar2] n {6 > 0}) ^ 0, 

where CAD denotes the symmetric difference of C and D. In fact, if (|4.13|) does not hold, we 
have that there exists x G i°°(R^) with x^{b>a} = 0, x 7^ and 

l[o„AO]n{b>o} X weakly* in L°°(R^). 
As a consequence it would be 



bl[n„An]dx^ / &l[o„Aa]n{6>o} c^^; -> / bxdx > 

which is against H4.12(l . Since 

|meas(r2„ (1 {b > 0}) - meas(fJ n {& > 0})| < meas([17„Ar2] n {& > 0}) 0, 

we deduce that (|4.11|1 holds. 

Let us assume now that (|4.11|) holds. Let us set ii„ :— mo„. Choosing u„ as a test in (|4.1U|I 

we deduce that m„ is bounded in L^'''(]R^). Up to a subsequence we have that there exist $ G 
LP(M^M2) ^ £ LliW^) such that 

(4.14) Vu„lo„^$ weakly in LP (R^M^) 
and 

(4.15) u„lo„ ip weakly in L^(R^). 
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By the convergence of f2„ to Q in the Hausdorff complementary topology, we have that $ = \7ip on 
n. In fact let € C^i^, M^) and let U he a regular subset of Q such that supp(*) CU CU cn. 
Then U C for n large and 

Vu„ • '5 da; = — / u„div'I' dx. 
Ju 

Let c„ be the average of m„ on [/. By Poincare inequality and Rellich Compact Embedding of 
W^'P{U) into L^{U) we get up to a further subsequence 

(4.16) K-c„)^u strongly in LP (f/). 
In particular the convergence is strong in L^{U) and 

(4.17) / ^-^dx^- I udiv^dx. 

Ju Ju 

By (|4.15|) and (|4.16|) we deduce that c„ = w„ — — c„) converges to some c € M. We conclude 
that u = if — c and by H4.17|l we get 

<i> • dT = - / ip div^f dx 

u Ju 



which means that $ = V(/3 on fi. 

Notice moreover that (|4.1HI and 14.15|l imply that 



(4.18) lim / Unln hdx = / tphdx. 

n-^+oo 7jj2 " Jn 

In fact, since h = bf + g and supp((7) C fi, it suffices to check that 

lim / u„lsi bf dx = / (pbf dx. 
Since ln„n{fc>o} ^ ^nn{b>o} strongly in ^^(R^) we have 



lim / u^ln^bfdx^ lim / w„lo„^/lo„n{b>o} c^a^ = / vbfl^r^{b>o}dx. 
Let us prove that ip — un on fl. In fact, for every v G £^'''($1), by monotonicity we have that 
(4.19) / [A{x,Vvln)-{S7vln-Vu^lnJ + B{x,vln)ivln-UnlnJ]dx 

> / [A(x, Vu„lo„) • (Vwlo - Vu„lo„) + S(a;,u„la„)(t;ln - u„lo„)] dx, 
Jr^ 

where 



We claim that there exists Vn G i^'''(ri,i) such that 



B{x,0 ■.= b{xmP-^^-hix). 

(4.20) Vt)„lo„ ^ Vvln strongly in LP (M^R2) 
and 

(4.21) i'„ln„-^'yln strongly in LP (M^). 
Notice that 

(4.22) lim / hvnln dx = hvl^dx. 
In fact, because of H4.21|l . we have 



n — > + oo 



lim / bfvnln„dx = / bfvl^dx. 



Moreover, if (7 is regular and with supp(.g) C [/ C J7 C we have by Poincare inequality 



g{vn - v) dx 



< II5|Ilp'(c/)I|Vu„ - Vu||i,p([/,R2) -> 0. 
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Since U C f2„ for n large enough, we conclude that (|4.22|) holds. 

Using Vn — u„ as a test function in H4.10|l we can rewrite the right-hand side of (|4.19|) as 



(4.23) / [A{x,yunlnJ-{yvln-yunlnJ + B{x,UnlnJ{vln-UnlnJ]dx 

= / [^(a;, Vw„ln„) ■ (Vuln - Vu„lf2„) + B(x, M„lf^„)(t;ln - w„ln„)] da;. 

Since A{x,Vunln„) is bounded in (M^,R^) and |un|^"^w„lf2„ is bounded in (K^), passing 
to the hmit in H4.19|l and in l|4.23(l . by claims H4.20(l and H4.21II . and in view of (|4"l8(l and of the 
fact that A{x, 0) = 0, we obtain 



n 
> 



(4.24) / [A{x, Vv) ■ ( Vw - V(y5) + B(x, v) {v - (p)] dx 

lim / [A{x, Vw„ln„) • (Vwlji - Vw„ln„) + B{x,Unln,J{vln - Wnlo„)] dx = 0. 

'^ + °° JR2 

Taking v = (p ± ez in (|4.24ll . with z G L^'^(f7) and e > 0, dividing by e, and passing to the limit 
as £ ^ 0, we obtain that ip = uo in Q. 

Let us now prove that (p — un^n, ^ ~ Vuolo, and that the convergences in (|4.14|l and H4.15|l 
are indeed strong. Let us take v :— un in H4.23(l . Since A{x, 0) = we obtain that 



lim_ / {A{x, Vu„lo„) - A{x, Vuolo)) ■ (Vu„lo„ - Vm^Io) dx 

+ / b{x){\un\P^'^Unln„ - \un\^~'^unln){unin^ - unlfi)]dx = 0. 



By monotonicity we get that each integral tends to zero. Now, the strong convergence of Vw„lo,^ 
to VuqIq in L^'(M^,R^) and of w„la„ to uqIsi in L]J(M^) is a consequence of ^1 Lemma 2.4]. 

In order to conclude the proof, we have to prove our claim on the existence of w„ G L^'^(S1„) 
satisfying H4.20|l and (|4.21|l . By Proposition 14.21 every function z e L^'^(ri) is strong limit of 
a sequence of functions z„ G -^b'^(^n) (with the usual extension to zero outside fl and fin re- 
spectively). Since fl G Ap{M.^), by Proposition 14.11 v is a strong limit in L^'^(ri) of functions in 
W^''^{fl), which in particular are in £^'^(0). Then the strong approximability for v that we need 
follows easily using a diagonal argument. □ 

4.2. Nonlinear optimal cutting problem. In this subsection, we apply Theorem 13.81 to the 
problem of optimal cutting for a membrane governed by a nonlinear energy. 

Let r2 C be open, bounded, and with a Lipschitz boundary, and let xi,X2 G fl. Let us set 

IC{fl) := {K C_ fl : K is compact and connected with xi,X2 G K}. 

Let / : R-^ X ^ [0, +oo] be a Caratheodory function such that f{x, 0) = and satisfying the 
following growth estimate 

(4.25) «|er</(2:,e)<a(ier + l), 

where a > and p g]1, +oo[. 

Let g G W^-p{R^). For every K G K-iTl) let us set 



(4.26) 



£{K) := inf \ [ f(x, Vu) dx : u e L^^P{n \ K),u = g on dn\ k \ , 
[Jn\K J 



where L^ P{il. \ K) is the Deny-Lions space defined in (|4.1|l with 6 = 0. Notice that it is natural 
to consider displacements in i^'^'(f2 \ K) because the energy involves only Vu, and so we cannot 
expect to control the L^-norm of u. Moreover, notice that the boundary condition on dil \K is 
well defined since Vl is Lipschitz and u G W^^P{fl fl Br{x)) for every x G dfl \ K and r such that 
Br{x)^K = 0. 
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We have that £{K) can be rewritten as 

(4.27) £{K) ^ mf \ / f{x, Vu) dx : u <E W^P{n \ K),u = g on dn\ K \ . 

\^Jn\K J 

This is due to the density of W^'P{Q \ K) in L^'P{^ \ K), but a Httle care should be paid 
for the boundary condition. In particular, denoting by Tm the truncation operator Tmu :— 
niin{niax{M, —M}, M}, we have for every u G L^'P{n \ K) 

I f{x, Vu) dx — lim / f{x, VT^u — X/TmS + V.g) dx 
Jn\K '^^^+°°Jn\K 

so that (|4.27|l holds. 

The optimal cutting problem consists in finding the "cut" K G ]C{Q) which maximizes £ among 
all admissible cuts, i.e., to solve the problem 

(4.28) max_ £{K). 

The existence of an optimal cutting has been established by Bucur, Buttazzo and Varchon in |S] 
in the case p = 2 and with a quadratic energy density /(x,^) = ■ ^. In view of Theorem 13.81 
the existence of an optimal cut can be proved also in the nonlinear case 1 < p < 2. The following 
result holds. 

Proposition 4.5. Let 1 < p < 2. Then problem H4.28|l has a solution. 

In order to prove Proposition 14.51 we need the following lemma which is based on Theorem l3.8l 
and on the first condition of Mosco convergence for Sobolev spaces W^'"^ proved in 0. 

Lemma 4.6. Let 1 < p < 2, and let {Kn)n(EN be a sequence in IC{n) converging in the Hausdorjf 
metric to K . Then for every u G W^'^{Q,\K) with u = g on dfl\K there exists Un G W^'P{fl\Kn) 
with Un — g on dfl \ Kn such that 

Vu„lo\K,. Vuln\K strongly in LP(M^,M^). 
Proof. Let B be an open ball containing f2. Let us consider 

{u in Q\K 
g inB\n. 

Notice that u G W^^p{B \ K). By applying TheoremEHto {B \ K) G Ap{^^) we have that there 
exists Vh G W^ '^{B \ K) such that 

ih^u strongly in T4^i'P(B\iC). 

By Proposition l4.2l fwith & = 0, so that the convergence of the measures is automatically satisfied), 
for each G N there exists -DJJ G W^''^{B \ Kn) such that 

Vw;J1b\x„ ^ ^ihlB\K strongly in L''{m.^,m?). 

Using a diagonal argument, and since the L^-convergence is stronger than L^-convergence on 
bounded domains {1 < p < 2), we can find Un G W^'^{B \ Kn) such that 

(4.29) \/UnlB\K„ -^'^u1b\k strongly in LP (R^R2). 

The functions u„ do not satisfy a-priori the required boundary condition, so that we have to 
modify them. We follow here an idea due to ChamboUe Up to adding a constant to t2„, we 
can assume that 

u„ ^ g strongly in W^'P{B \ U). 
Let us consider gn := {g — Un)\B\n- have that gn G W^'P{B \ fl) with g„ strongly in 
W^'P{B \ fi). Let E denote a linear extension operator from W^ p{B \ fl) to W^ p{B): such an E 
exists because has a Lipschitz boundary. Let us set 

Un ■= (Un + Egn)\n\K„- 
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We have m„ G W^'P{^1 \ K^) with 7i„ — g on Sfi \ Kn. Moreover, since Egn strongly in 
W^'P{B), in view of H4.29|l . we deduce that (u„)„eN is the required sequence. □ 

We are now in a position to prove Proposition ^3] 

Proof of Proposition \4-5[ Let {Kn)neN be a maximizing sequence for the optimal cutting problem, 
i.e., 

£{Kr^) snp{£{K) : K e /C(n)}. 

Up to a subsequence, we can assume that Kn K in the Hausdorff metric. We have that K is 
an admissible cut, that is K G IC{n). 

By Lemma Ol for every u G M^i'P(f7\i^) with u = g on dn\K there exists m„ G M^i'P(rj\i^„) 
with Un = g on dQ \ Kn such that 

Vu„lo\J^„ -> VMlo\/f strongly in LP{R'',R^). 

Since f{x, 0) = 0, we deduce 

/ f{x,\/u)dx— lim / f{x,Vun)dx> lim £{Kn)- 

Jn\K n->+oo n^+oo 

Taking the infimum over all admissible u we get 

£{K) > lim £{Kn) 

n — > + oo 

SO that K is an optimal cut, and the proof is concluded. □ 

Remark 4.7. In the proof of Proposition l4.5l it is not clear if meas(i^„) meas(if ) for n — > +oo. 
As a consequence, the result of Dal Maso, Ebobisse and Ponsiglione ^H] cannot be applied to 
recover the approximability of gradients of functions in W^'P{^\K) through gradients of functions 
in W^'P{^1 \ Kn) that we need. It seems essential to use the approximability of the gradients for 
the relative VF^'^-spaces established in |5] and the density result given by Theorem 13.81 

Let us assume that f{x,^) is strictly convex in ^ so that the problem 

(4.30) min \ / f{x, Vu) dx : u e L^'P{n \ K),u ^ g on dn\ K \ 

[Jn\K J 

admits a unique solution Ufi\x G L^'P{fl \ K). In particular, in view of 14.26() we have 

^{K)=j f{x,Vun\K)dx. 
Jn\K 



The associated Euler-Lagrange equation is 
(4.31) 



divd^fix, Vuii\K) =0 inn\K 
d^f{x, Vufi\K) ■ 1^ = on K 
un\K ^ g ondn\K. 

We deduce that the following stability result for the Neumann-Dirichlet problem (|4.31|l holds. 



Proposition 4.8. Let 1 < p < 2, let K be a solution of (|4.28|l . and let (A'„)„gN be a sequence 
in IC{fl) converging in the Hausdorff metric to K . Then we have that Q.\ K is stable for H4.31|l 
along the sequence (fl \ Kn)nGN, that is 

(4.32) yun\K„ ln\A'„ ^ Vun\Ala\if strongly m LP{R^, R^). 

Proof. Choosing g as an admissible displacement in (|4.30|) . we get that VuQ\A'„lsi\i<'„ is bounded 
in LP(M^,M^) so that up to a subsequence we have 

^un\Kjn\K^ ^ ^ weakly in LP{R^,R^). 
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Since Kn K in the HausdorfF metric, there exists u e L^'P{il \ K) with u = 5 on dVL \ K such 
that <i> — Vu on ^l\K . Moreover by lower semicontinuity we have that 



(4.33) / f{x,\/u)dx< [ /(x,$)(ia; < liminf / f{x,Vun\Kjn\Kjdx 

= liniinf/ f{x,\/un\K„)dx. 
"^+°° Jn\K„ 

Let V e L^^P{fl \ K) with v ^ g oti dVl\K. By Lemma lOl there exists u„ £ L^'PiVl \ Kn) with 
Vn = .9 on 9il \ Kn such that 

Vw„lo\A'„ ^ Vvla\K strongly in LP(R2,R2). 

Then by we get 



(4.34) / f{x,\7v)dx^ lim / f{x,Vvnln\Kjdx= lim / f{x,Vvn)dx 
Jn\K "^+°oJr2 "^+°°Jo\x„ 



> liminf / f{x,\/uQ\j( ) dx > / f{x,Vu)dx 

"^+~Jo\A'„ " Jo\A 

SO that w = u^2\K- Taking v = u in (I4.34|) we obtain 



lim / f{x,\7un\Kjdx ^ f{x,Vua\K)dx. 
^^+'^Jn\K„ Jn\K 

By 0] we conclude that H4.32|l holds, and the proof is concluded. □ 

5. Appendix: the density result for the symmetrized gradients 

The aim of this appendix is to show how our approach to density explained in Section U] can be 
employed to prove the density of the Sobolev space W^''^ in the spaces LD^'^ of two dimensional 
elasticity. The case p = 2 is the really interesting one, and has been proved by ChamboUe in |13|: 
using this density, he proves existence for the Cantilever Problem and for the evolution of brittle 
fractures in the context of planar linearized elasticity. Our approach provides a different proof of 
ChamboUe's result, and covers also the case 1 < p < 2. 

In order to make the context precise, let O be an open subset of M^. For I < p < +00, let us 

set 

LD^'PiQ) {u e W.'-Pin^R^) : eiu) e LP(f7, M.^^^)} , 

where e(u) := (Vu + (Vu)"^)/2 denotes the symmetrized gradient of u, and M^^j^ denotes the 
space of 2 X 2 symmetric matrices. Clearly W^'P{^},M.'^) C LD^ P{n). If H. is Lipschitz, by means 
of Korn's inequality, it turns out that LD^'P{U) coincides with W^'P{il.,M.'^), while if is irregular, 
the inclusion can be strict. 

The main result of this section is the following. 

Theorem 5.1. Let 1 < p < 2, and let 17 C be a bounded open set such that ff^ has a finite 
number of connected components. Then for every u G LD^'P{fl) there exists u„ G W^''^{fl,M.'^) 
such that 

e{un) e{u) strongly in L^(il, M^^j^). 

Proof. Let Ki, i — 0,1, ... ,m, he the connected components of il^, where Kq is the unbounded 
one. Let us consider the space 

H {e{v) : v £ H\n,R')} C LP{n,Ml^^J, 

where on M^yj^ we consider the scalar product A : B := tr{AB'^) — J^i j o-ijbij- 

In order to prove the theorem, it suffices to check that for every u G LD^'P{^1) we have 

e{u) e Tl, 

where the closure is taken in the L^-norm. 
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We employ a functional analysis argument, namely that H = (H^)-^, where (•)"'" denotes the 
orthogonal in the sense of Banach spaces. So our strategy is the following. Firstly we characterize 
H^, and then we check that e{u) is orthogonal to every element of 

Step 1: Characterization of H^. We claim that 

(5.1) i/-^ |Hess((p) : (p e W^'^' (R^), tp is linear on Ki for i = 0,1, ... ,mj , 

where Hess(iy9) is defined as 

By linearity of ip on Ki we mean that there exist c,; G and 6,; G M such that (notice that 
Wo'P'(R2) c CHK')) 

(5.3) </5(a;) = Ci ■ x + bi for x G Ki. 

Since (p G W^q (R^), we clearly have cq = and bo ~ 0. 

Let us check H5.1(l . Let * G i^'(f^, M^yj^) be an element of H^, where p' := p/{p — 1) is the 
conjugate exponent to p, with 

viz = f ^1 ^2 

V V'2 V'3 

This means that for every v G i?^(ri,R^) we have 



n 



m 

From Lemma [3.41 we conclude that there exist 0i, 02 G W^-'^'^ (R'^) and G M^, i = 0, 1, . . . ,to 
such that 

V01 =i?(V'l,V2), V(/)2 =-R(V'2,V'3), 

(5.4) (</>!, 02) = c, on /C, for 1 <p<2, 
and 

(5.5) (01, 02) = Ci C2-q.e. on i^i for p = 2, 

where R{a, b) :— {—b, a) denotes a rotation of 90 degrees counterclockwise. We can assume co ~ 0, 
hence 0i,02 G Wq'"' (R^). 
Let us set 

(5.6) (01,02) G<''''(RMR2). 

Let be a disk centered at the origin and such that fl C D. For every v G H^{D,M.'^) we have 
that 

/ $-Vwdx = - / {div^)vdx = - / (9i0i + 9202)^ rfa; = - / {--tjj2 + 'ip2)v dx ^ 0. 
Jd Jd Jd Jd 

Using again Lemma [3.41 we get that there exists ip G W^'^ (M^) with ip = on R.'^ \ D, and such 
that 

V(p = i?$= (-02,01 ). 



^Jf : e{v) dx = 0. 

Choosing w G of the form v := (ui,0) with vi G H^{n), we deduce that for every vi G H^{n) 

/ (^^1^2) • Vwi da; = 0. 
Jn 

Similarly we deduce that for every V2 G H^{il) 

(-02,V'3) • Vw2 da; 0. 
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In view of (|5.6I) . we conclude that £ W^o (M^). Since > 2, by Sobolev Embedding Theorem 
we have that G C^(R^), so that, by Lemma [2. II from (|5.4|l and (|5.5(l we get (up to replacing a 
with i?Ci) 

(5.7) V<^ = on iCi. 

By construction we have that \1/ = Hess((y9). In order to complete the proof of claim (|5.1|l . we 
need to check l|5.3|l. Let us consider 

ipi{x) := ip{x) ~ a - X. 

By (15.7(1 . we clearly have that Vi^i = on i^i, i.e., Ki <Z d, where Ci is the set of critical points 
of ipi. By Sard's Lemma we have that 

meas(v3i(C,;)) = 0. 

Since ipi{Ki) is connected, and u\eas{ipi{Ki)) ~ 0, we conclude that ipi{Ki) ~ {bi} for a suitable 
bi £ M, so that 1(5.3(1 is proved. 

Step 2: Checking the orthogonality condition. Let u G LD^ ^(rj), and let ^ G iJ^. We 
have to check that 

* : e(u) = 0. 

According to ((5.1(1 . let ^ = Hess((/?) with ip G Wq'^ (M^) satisfying ((5.3() . Let us consider G 
C°°(R^) and G C^iM."^), ? = 1, 2, . . . , to such that = 1 on a neighborhood of Kq, f ^ = 1 on a 
neighborhood of Ki, and 

supp(^h) n supp(^fe) = ior k. 
By [3 Theorem 9.1.3] we can find ipl G C°°(K2) with 

(Pn{x) — Ci ■ X + bi on a neighborhood of Ki, 

and such that 

Strongly in VF^''''(R2). 

Let us set 

(m \ rn 

Clearly we have that 

(5.8) 'Pn-"f strongly in VF2'P'(R2), 
and 

(5.9) Hess(<p„) = on a neighborhood A„ of fi^. 

We can assume that \ An is regular. Then by means of Korn's inequality we have that 

(5.10) u&W^^P{n\A^,R^). 
By ((5.8() and ((5.9() we conclude that 

vl^ : e{u)dx— / IIess(93) : e{u) dx = lim / Hess((y9„) : e{u)dx 



= lim / IIess(iy9„) : e{u)dx. 

By ((5.10(1 and since Hess(iy9„) is symmetric, we deduce that 

(5.11) / ^ : e{u)dx— lim / Hess(iy9„) : Vudx. 
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Notice that the rows of Hess(iy9„) are divergence free in il \ and with nuU trace on 9(17 \ An). 
Integrating by parts in (|5.11|) . we get 

* : e{u) dx = 0, 

so that the proof is concluded. □ 

Remark 5.2. In his proof of the density of W^''^{fl) in LD^-^{n), ChamboUe ^21 considers 
LD^'^{il) (up to functions u such that e{u) = 0) as a Hilbert space with scalar product {u,v) := 
e{u) : e{v) dx, and proves that 

{e{u) : u e W^^^i^)}^ = 0, 

where (•)^ is the orthogonal in the sense of Hilbert spaces. In this framework, the function 
appearing in our Step 1 is of the form 5* = e{v) for some v G LD^'^{U), and the same analysis 
implies that e{v) = Hess((/3). As a consequence we get A^(/j = (93 is usually called the Airy 
function). ChamboUe uses some PDE and capacity arguments to show that = Q \n the case J7 
is simply connected, and then proves the general case by reduction to the simply connected one. 

In our case, we cannot employ PDE arguments, because we consider LD^'P{U) as a natural 
subspace of L^i^, M^y,^), and this seems unavoidable in the case 1 <p <2. As a consequence our 
function ip does not satisfy A^(^ = 0, and we must work out an approximation of as in Step 2. 

Remark 5.3. In order to follow the arguments of Step 2, it suffices to approximate ^I' — IIess((y9) G 
by '^n G (17, l^sym) which are null on a neighborhood of 17^ and whose rows are divergence 
free. This is obtained taking '■= Hess((y5„), where G Wq'^ (M?) is such that 

ipn~^^ Strongly in (R^), 

and with linear on a neighborhood of 12^. This last constraint cannot be treated using ideas 
similar to Lemma 13.31 so that we used partition of unity (which requires Vl'^ with a finite number 
of connected components) and the approximation result [3 Theorem 9.1.3] (which requires 1 < 
p < +00). 
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